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ABSTRACT 

We show how analysis of a quasar high-magnification microlensing event may be used to construct a 
map of the frequency-dependent surface brightness of the quasar accretion disk. The same procedure also 
allows determination of the disk inclination angle, the black hole mass (modulo the caustic velocity), 
and possibly the black hole spin. This method depends on the validity of one assumption: that the 
optical and ultraviolet continuum of the quasar is produced on the surface of an azimuthally symmetric, 
flat equatorial disk, whose gas follows prograde circular orbits in a Kerr spacetime (and plunges inside 
the marginally stable orbit). Given this assumption, we advocate using a variant of first-order linear 
regularization to invert multi-frequency microlensing lightcurves to obtain the disk surface brightness as 
a function of radius and frequency. The other parameters can be found by minimizing in ^ fashion 
consistent with the regularized solution for the surface brightness. 

We present simulations for a disk model appropriate to the Einstein Cross quasar, an object uniquely 
well-suited to this approach. These simulations confirm that the surface brightness can be reconstructed 
quite well near its peak, and that there are no systematic errors in determining the other model param- 
eters. We also discuss the observational requirements for successful implementation of this technique. 

Subject headings: accretion, accretion disks — gravitational lensing — quasars: individual (2237-1-0305) 
— relativity 



1. INTRODUCTION 

Due to their great distance and small intrinsic size, it is 
not possible to obtain a resolved optical image of a quasar 
with current technology. The angular size of a quasar op- 
tical emission region is of order 10~^ , a scale so small as 
to require a baseline of several thousand kilometers to re- 
solve. Thus, until optical VLBI becomes practical, quasar 
structure will need to be probed by other, indirect, means. 
Reverberation mapping provides an instructive example 
of the difficulties of such indirect approaches, for it has 
proven difficult both to implement and to interpret. We 
believe that the subject of this paper, microlensing by stars 
in an intervening galaxy, is a more promising method. 

The circumstantial evidence that black holes power 
quasars is convincing: accretion onto black holes can be 
very efficient in converting rest mass energy to photons (up 
to 40%) or to bulk momentum (forming radio jets/lobes) 
in a compact region which has an effective temperature 
near where the quasar spectrum peaks. There is statis- 
tical evidence that quiescent black holes in the nuclei of 
galaxies at low redshift could be the remnants of quasars. 
And, quasars have properties very similar to Seyfert galax- 
ies, for some of which there is good spectroscopic evidence 
for a central black hole. The nature of the accretion flow is 
quite uncertain, however, although it is probably geomet- 
rically thin since angular momentum will support the ac- 
cretion flow against collapse and geometrically thick disks 
are by nature inefficient (if quasar disks were thick, there 
would be problems producing the huge luminosities ob- 
served within a reasonable mass budget). 

Attempts to constrain the character of the innermost 
accretion flow by means of spectral modeling have made 
little progress for several reasons. There are major sys- 
tematic uncertainties about fundamental issues (e.g., the 



vertical distribution of energy dissipation, the physics to 
include in radiation transfer solutions). In addition, any 
particular model depends on a sizable number of free pa- 
rameters (mass of the central black hole, accretion rate, 
viscosity parameter, inclination angle), so that parameter 
estimation is tricky. 

Analysis of a microlensing event is potentially a more 
powerful tool. Rather than guess a specific model for the 
accretion disk, the history of magnification in the event 
can be used to directly infer the disk surface brightness 
as a function of radius and frequency. The only assump- 
tions required are that the continuum emission surface is 
geometrically flat, the material forming it follows circular 
orbits, and general relativity determines the dynamics of 
both the matter and the photons. 

To apply this technique requires study of a particular 
gravitationally-lensed quasar. Many are now known, and 
some appear to undergo occasional fluctuations due to in- 
dividual stars in the lens galaxy magnifying the quasar, 
a phenomenon referred to as microlensing. The Einstein 
Cross is particularly well-suited to this problem for a num- 
ber of reasons that we will detail in §1.3. Several authors 
have attempted to constrain the character of the accretion 
disk in this system by comparing predicted lightcurves to 
the data compiled during microlensing events (Jaroszyhski 
et al. 1992, Ranch & Blandford 1991, Jaroszyhski & Marck 
1994, Czerny et al. 1994). However, the results have all 
been somewhat inconclusive due to the ordinary spectral 
modeling difficulties described above. 

The remainder of this paper is organized as follows: In 
§1, we describe how a caustic crossing can be used to per- 
form the mapping, and discuss why the Einstein Cross is 
such a suitable target for this sort of study. In §2, we 
discuss the model in detail, considering our assumptions. 
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the inversion technique, and error propagation. In §3, we 
present a variety of simulations of caustic-crossing events, 
and demonstrate that we can measure most model param- 
eters. We show some examples of the accuracy with which 
the intensity at the accretion disk may be determined, and 
discuss how the reliability of the results depends on the 
quality of the observational data, on the character of the 
monitoring, and, to some degree, on the character of the 
regularization scheme. In the final section, we discuss the 
results and present our conclusions. 

1.1. Smooth thin disk 

For a standard thin accretion disk with constant accre- 
tion rate and no advcction of heat, the energy generation 
per unit area as a function of radius is given by 



i/=10i» Hz 



3 GMbhM 



RR{r). 



(1) 



where Mbh is the mass of the black hole, M is the ac- 
cretion rate, r is the radius within the disk, and Rj^ is 
a correction factor that combines outward advection of 
energy associated with the angular momentum flux and 
relativistic effects (Page & Thorne 1974, notation from 
Krolik 1998). Rjf is a function of the black hole spin 
Us = a/MsH- 

Though equation describes the functional depen- 
dence of the energy released with radius, it does not spec- 
ify whether this energy is thermal or mechanical, and does 
not specify the dissipation as a function of height within 
the accretion disk. The appearance of a standard thin ac- 
cretion disk can vary significantly depending on how and 
where the energy is released. To further our understand- 
ing, it would be very desirable to actually measure the 
local spectrum. 

One potential method to achieve this is to observe a 
quasar during the sort of high-magnification microlens- 
ing event that occurs when a caustic crosses the source 
(Grieger et al., 1988 and 1991, Gould & Gaudi 1997). 
Grieger et al. (1991) advocate inverting the microlensing 
lightcurve to obtain the one-dimensional surface brightness 
of the quasar, Pv{x) — J Iu{x,y)dy, where Iv(x,y) is the 
specific intensity of the quasar at sky coordinates {x,y). 
Their method is quite elegant, but relies on first-order 
regularization (the assumption that the quasar profile is 
smooth). This assumption is problematic for black hole 
models since the energy release increases rapidly towards 
smaller radii, and relativistic effects can cause sharp peaks 
in the quasar profile. These sharp features are smoothed 
out when this method is used (see Figure 1). 

We believe a better approach is to solve instead for the 
surface brightness at the accretion disk, with the relatively 
benign assumptions that the disk is planar, axisymmetric, 
and isotropically emitting (in the fluid frame), and most 
importantly, likely to vary smoothly with radius. This ap- 
proach provides a better inversion of the quasar profile if 
the disk assumptions are correct. This method recovers 
the spectrum as a function of radius at the accretion disk, 
which can then be compared with disk atmosphere mod- 
els or other spectral modeling. In addition, it can help 
constrain disk parameters, such as the inclination angle, 
without relying on a specific accretion disk spectral model. 
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Fig. 1. Plot of recovered one-dimensional disk profile as 
a function of impact parameter for a face-on accretion disk 
using the Grieger et al. (1991) technique. Solid line is the 
model profile; dotted line is recovered profile. The param- 
eter a is distance across the source perpendicular to the 
caustic line in units of gravitational radii. 

1.2. High magnification microlensing event 

The spectrum of a quasar as a function of time is a 
convolution of the lensing magnification with the surface 
brightness of the quasar. Near a fold caustic, two bright 
images merge (Schneider et al. 1992). The resulting mag- 
nification has the specific form: 



A(x,t) = Ao + 



K 



\Jx- Vcit - to) 



Q[x- vc{t - to)] , (2) 



where Aq is the magnification due to the additional im- 
ages, X is the position on the quasar plane perpendicular 
to the caustic, K is the strength of the caustic [units of 
(distance) ^/^], Vc is the speed of the caustic (assumed to 
be positive), Q is the step function [8(a;) = 1 for a; > 0, 
otherwise] , and to is the time at which the caustic crosses 
the X — point measured relative to the center of the 
lightcurve (Schneider & WeiB 1986). This equation is valid 
whenever the source size is small compared to the Einstein 
ring radius of the microlensing star. If that ratio is small, 
there is only a small probability of any of a number of 
problems that might invalidate equation [||. The list of 
potential problems includes: 1) the possibility that the 
source is projected behind a cusp (where the caustic curve 
discontinuously changes direction) or behind the crossing 
of multiple fold caustics; 2) significant curvature in the 
caustic; 3) variation of K along the caustic on the scale 
of the source. Grieger et al. (1988) demonstrate that a 
source of size smaller than ~ 10% of the Einstein radius 
of the typical lensing mass is necessary for the second and 
third assumptions to be valid. 
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1.3. Microlensing Laboratory: The Einstein Cross 

The Einstein Cross is a quasar (z=1.695) lensed into 
four images by a nearby (z=0.0394) barred spiral galaxy 
(Huchra et al. 1985). The luminosity distance to the 
quasar is Dl = 3 x lO^^h^^ cm, where /175 = Ho/{75 
km/s/Mpc), assuming = 1, A = 0. 

This quasar has been observed to undergo fluctuations 
due to microlensing roughly once per year (Irwin et al. 
1989; R. Webster, private communication). It is partic- 
ularly well suited for studying microlensing because the 
lensing galaxy is nearby. This happy coincidence makes 
the time delays between the four images all less than a 
day, so that it is easy to distinguish intrinsic variability 
from microlensing variability. It also makes the stellar ve- 
locities projected onto the source plane quite high, so that 
microlensing is frequent, and also makes individual events 
comparatively brief (only about a month, in contrast to the 
years to decade timescales characteristic of more distant 
lenses). In addition, the Einstein radius projected onto the 
quasar plane is quite large compared to the quasar size, 
validating the assumptions made in the previous subsec- 
tion, and also making the microlensing variations strong. 

The models of Witt et al. (1993) suggest that for image 

A, 

<^)-K 5.8xToO ^^^-"'' 

where = 5.8 x 10^^cm(m/0. 2)^/^/1^5^^^ is the Einstein 
radius projected to the quasar plane (Schneider et al. 
1992), TO is the typical mass (in Solar units) of a star caus- 
ing microlensing, and Mg — AIbh /W^Mq. In equation 
(3), and in the rest of the paper, we adopt rg = GM/c^ as 
the unit of distance. For concreteness, we will use K = ^va 
our simulations. The parameter Aq can be approximated 
by {Aq) — 1(1 — cr)^ — 7^1"-^ where 7 is the shear (Witt 
et al. 1993). For image A, the estimated range of the mi- 
crolensing parameters is cr = 0.3 — 0.4 and 7 = 0.4 — 0.5, 
giving Aq = 3 — 9; we use Aq — & in our simulations. 

2. MICROLENSING OF AN ACCRETION DISK 

2.1. Caustic crossing 

As the caustic crosses the quasar, the observed 
lightcurve is a convolution of the magnification, equation 
P] and P^{x): 

/oo 
A{x,t)P^{x)dx. (4) 
-00 

Note that if the flux within a waveband can be measured 
outside of the caustic and subtracted off, then the depen- 
dence on Ao disappears, and the K parameter becomes de- 
generate with an arbitrary scaling of P^ix). As discussed 
by Grieger et al. (1991), equation Q can be inverted us- 
ing regularization to flnd P,j(x). Similar techniques have 
been used for measuring the limb-darkening of stars during 
galactic microlensing events (Gaudi & Gould 1998, Albrow 
et al. 1998). 

2.2. Black Hole geometry 

Near a black hole, relativistic effects cause Doppler 
beaming of the emitted radiation, gravitational red shifts, 



and bending of photon trajectories. To image the surface 
of an accretion disk, these relativistic effects must be ac- 
counted for using a relativistic transfer function (defined in 
Cunningham 1975). To compute the transfer function (see 
equation Q below), we make several simplifying assump- 
tions: (1) the accretion disk is thin, i.e. h ^ r; (2) the 
gas follows prograde circular orbits outside the marginally 
stable radius rms, and undergoes freefall within r^s with 
constant angular momentum and energy equal to those ob- 
taining at rms ; (3) the disk is fiat and lies in the equatorial 
plane of the black hole; (4) the gas emits isotropically in its 
rest frame, i.e., there is no limb darkening in the accretion 
disk atmosphere. The first two assumptions are appropri- 
ate if pressure gradients cause forces much smaller than 
the gravitational force in the z and r directions, respec- 
tively. This condition is not met in advection-dominated 
accretion flows or slim accretion disks (Beloborodov 1998). 
In the case of slim accretion disks, the orbital frequency 
is nearly Keplerian, and deviates by less than 20% when 
Mc^ / Lsdd < 1000; however, the disk scale height can be- 
come a large fraction of the radius, which changes the emit- 
ted angle of radiation relative to the disk normal and can 
cause shadowing which we do not take into account. The 
third assumption is inappropriate if the disk is warped; 
however, Bardeen-Petterson precession (1975) can align 
the disk and black hole by the time the gas reaches the 
inner radii. The fourth assumption is a simplification for 
greater ease in the inversion computation since the disk 
can be viewed from only one angle and thus at most one 
emitted angle can be observed at each radius/azimuth of 
the disk. For each radius, there is a limited range of emit- 
ted angles which are observed, so our inversion will give 
some sort of average of the intensities within that range. 
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Fig. 2. The upper curve shows the maximum fi^. — cos0, 
where 9 is the angle between the disk normal and the di- 
rection of the photon in the fiuid rest frame. The lower 
curve shows the minimum /Ze- The disk parameters are 
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i = 30° and a,, = 0.998. 

Figure 2 shows the range of emitted angles (in the fluid 
rest frame, /ie is the cosine of the normal to the disk) 
for a disk inclined at 30°. For a face-on disk, only one 
emitted angle is seen at each radius for all azimuths, so 
this assumption simply corresponds to mapping the spe- 
cific intensity of the disk at Heir). We could have assumed 
some limb-darkening law, but this is not warranted by the 
crudeness of the inversion technique. 

Figure 3 shows the disk geometry. The inclination angle 
of the accretion disk, = cosz, is z = 0° when the disk 
is face-on, and 90° when the disk is edge-on. The caustic 
crossing angle (j)c is measured with respect to the (a,/?) 
coordinates, which are defined so that the (3 coordinate 
lies parallel to the projection of the disk axis onto the sky 
plane (in Figure 3, the disk spin axis is pointing out of the 
page), with the black hole at the origin. We will use units 
of rg for the (a, (3) coordinates. 




Fig. 3. Geometry of the accretion disk. The disk axis 
points up out of the page, while the /3 axis lies in the page. 
The line parallel to the y axis is the caustic. 

The rotation of the accretion disk causes beaming of the 
radiation, leading to a hot spot on the approaching side, 
and a cold spot on the receding side. For an exactly face- 
on disk, the disk is symmetric, so no hot /cold spots exist, 
but a dip occurs inside the inner edge of the disk. In Fig- 
ure 4 we show {x) for a blackbody accretion disk at two 
frequencies. As the disk becomes more edge-on, the profile 
becomes more asymmetric as the Doppler aberration be- 
comes stronger. If the temperature of the disk decreases 
outwards, then the size of the hot spot will increase for 
smaller frequencies, and become less asymmetric, as can 
be seen by comparing Figure 4(a) and 4(b). Figures 4(c) 
and 4(d) show Pv{x) for different caustic crossing angles. 



showing the hotspot is oblong. 
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Fig. 1. Plots of the disk 1-D profile for a disk with 
Mg = 1, rn = 1, and = 0.998. Figures (a) and (b) 
compare inclination angle: i = 0° (solid), 30° (dotted) 
and 60° (dashed) for i/ic = 0. Figures (c) and (d) compare 
(j)c = (sohd), 30° (dotted), 60° (short dashed), 90° (long 
dashed) for i = 30°. The frequencies are in the quasar rest 
frame. 

2.3. Prediction of the Lightcurve 

With these assumptions and parameter choices, the ob- 
served flux may be predicted from the run of intensity with 
radius: 

F{i^o, t) = j dadpA [x{a, 0),t] g^I{uo/g, re), (5) 

where g = Vo/ve is the redshift between the observer 
and the emitter and I{ue,r) is the specific intensity at 
the accretion disk (we have assumed it is independent 
of the emitted angle). The relation between the caustic 
coordinate and the black hole coordinates is x(a,/3) = 
a cos 0c + Psin^c (see Figure 3). The physical variables 
Pe and re may be related via a Jacobian to a and /3 through 
their functional dependence on redshift g{iXo,<is), the in- 
clination angle of the disk /io, and the black hole spin a^. 
The magnification A depends onvc, to, 4>ci K, and Aq, for 
a total of seven model parameters. 

We compute the transfer function by shooting rays from 
infinity at a grid in (a, (3) until they cross the equatorial 
plane. The computational method is based on Rauch & 
Blandford (1994), and is described in Agol (1997). We use 
a nested grid of rays that is more finely sampled towards 
the center to resolve the inner parts of the accretion disk 
in greater detail. For a given observed frequency Vq, we 
compute Vf. as well as the emitted radius at each (a,/3), 
and interpolate these on the (pre-specified) grid of radii 
and frequencies at the accretion disk. 



Agol & Krolik 



To compute a normal transfer function, we would then 
simply sum over the grid. In this case, we multiply each 
ray by a further factor that describes the magnification 
due to the microlensing at any particular time, and then 
sum over the grid. This procedure may be summarized 



ih^,? (7) 



F,, 



T T 



(6) 



where i = I, Np {Np = NtX = ^ observed frequencies 
X # observed times) labels each measured observed fre- 
quency/time, and j — l,Ni [Nj = ^ emitted frequencies 
X # emitted radii) labels each emitted frequency/radius 
pair. The matrix T contains the integration and interpo- 
lation factors. 

2.4. Regularized Inversion 

Attempting to directly invert equation |6| for I given an 
observed set of F is impossible since the matrix T is gen- 
erally singular, so that noise in the lightcurve is magnified 
strongly during inversion. This fact requires the introduc- 
tion of some sort of a priori knowledge in order to make 
inversion feasible. Regularization is a particularly useful 
way to do this, as discussed in Press et al. (1992), because 
the "prejudice" injected into the solution is usually rela- 
tively benign and also relatively controllable. The essence 
of the linear regularization method is to minimize both the 
deviation of the model from the data, and also the devia- 
tion of the model from "smoothness," as defined by some 
sort of differencing operator. In our case, at any given fre- 
quency, we expect the emitted intensity I^, to be smooth 
as a function of radius, but not the one-dimensional profile 
P.- 

In our specific implementation of the method, we also 
impose several other restrictions on the solution. We ex- 
pect that the emitted intensity in the fluid frame dimin- 
ishes as the black hole event horizon is approached; we 
therefore require Ii^ir) to approach zero as r approaches 
Tg. In fact, because those regions are so strongly red- 
shifted from almost any inclination angle, the intensity 
in the fluid frame is almost completely unconstrained by 
the data, so physical assumptions have a very strong im- 
pact on the solution in this region. Similarly, we also re- 
quire 1 1, to approach zero at very large radii, for there is 
little energy available there to dissipate. Particularly at 
low frequency, it may sometimes be desirable to relax this 
constraint. In addition, we would like the inverted inten- 
sities to be positive definite; to achieve this, we maximize 
Aj — i^ihwij/Ni, where Wi is a weighting factor. The 
most appropriate weighting is Wi = rf since the radiat- 
ing area associated with each logarithmic radius interval 
scales as rf. Finally, for any choice of grid, there will al- 
ways be some radius/emitted frequency pairs that are not 
constrained by the data because Doppler shifts push Vq 
outside the observed region (these are the intensities for 
which the corresponding column of T is all zeros). In or- 
der to prevent those frequencies from contributing to the 
smoothing condition, we require that the associated inten- 
sities be zero. 

Combining all these considerations leads to the following 
regularization operator: 



B 




b:ifT,j=OV»} 

where is the number of frequency grid points at the 
accretion disk. The first sum describes the smoothness as 
a function of radius at each frequency, the two isolated 
terms give the boundary conditions at the innermost and 
outermost radii, and the last sum is the factor encourag- 
ing minimization of those Ij unconstrained by data. Note 
that the specific intensity vector is ordered with all the 
frequencies at one radius grouped together, so that li and 
li-N^ give the intensity at the same frequency, but ad- 
jacent radii. The smoothing operator B has the useful 
property of providing a model-independent measure of the 
"smoothness" of different solutions, due to the normaliza- 
tion by Al . Other differencing schemes might also be used; 
in the examples we have explored, it makes little difference 
to the outcome. 

Although the regularization condition is designed to be 
relatively innocuous, no such injection of prejudice can be 
altogether free from consequences (we will discuss the ef- 
fect of our particular choice in §3.2). We stress that the 
details of the regularization condition are always subject 
to "tuning" in the light of either theoretical expectations, 
or, better, the implications of real data. 

To solve for I, we minimize the following function: 



/(I) 



1 



x' + XB 



with 




(8) 



(9) 



where A is a constant, and ai is the error on Fi. We start 
with a direct solution of df/dl = 0, setting Aj ~ 1, 

I = ^M^M + -^H^ M^G = QM^G, (10) 

where My = Tij/a-j, Gi — Fi/cji, and H is defined such 
that B = 1 m 1. We then update Aj from the solution, 
and iterate until Aj converges. In some instances, Ai can 
become negative after an iteration. If this happens, we re- 
set Ai to be 0.1 times its value at the previous iteration, 
and recalculate the step. 

We assume that the magnification outside the caustic, 
Aq, can be measured from the lightcurve, and subtracted 
off. Then, the parameter K (the caustic magnification fac- 
tor) is completely degenerate with the disk surface bright- 
ness since a decrease in magnification corresponds to an 
increase in the surface brightness of the source. Conse- 
quently, we can determine the shape of the surface bright- 
ness profile, but not its absolute level. 

Several other parameters also remain to be determined 
after the direct inversion for the surface brightness profile. 
We call them collectively C, = (to, fc, 0c, ^s)- To find 
them, we fix A and compute over a coarse grid in this 
five-dimensional parameter space. Starting from the C, giv- 
ing the smallest i^i this grid, we refine our estimate of 
these parameters using the Levenberg-Marquardt method. 
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We compute the partial derivatives of with respect to 
the axes in space by finite differences. If a parameter 
with boundaries goes out of bounds, we fix it at the value 
where it went out of bounds, and keep it fixed throughout 
the rest of the minimization. This generally occurred with 
Us when it was near or 1 , and for ^ and 4>c when the disk 
was face-on. Fixing this improved estimate for the best- 
fit ^, we increase A and re-solve for the surface brightness 
profile until = ^f- 

The number of degrees of freedom against which to com- 
pare is not clearly defined for several reasons. One is 
that many of the model parameters are not entirely free; 
for several (^, Mbh, Vc) there are prejudices or constraints 
from other experiments. Another reason is the variable 
weight given the smoothing constraint vis-a-vis the data, 
as we are minimizing -\- XB rather than x^- In the 
limit of large A, there is effectively only one free parame- 
ter for each frequency in the fit to the 7^; in the limit of 
A = 0, there are as many free parameters as there are grid 
points. Given these considerations, Np is an upper bound 
to the true number of degrees of freedom; by raising A 
until x^ = ^F, we ensure that we do not overfit the data. 

2.5. Errors 

The word "error" has several different meanings in this 
context, and it is important to distinguish them. First 
of all, the errors in the inferred intensities have different 
properties from the errors in the model parameters. Sec- 
ond, both are potentially subject to systematic error as 
well as random error. 

We will begin by estimating the random error in the 
intensities I. Formally, we may say that 



6If 



E 



dFi 



(11) 



where 



dh 



J2kQikTk]/a-'^ (see equation [|o|), and 



assume the fluxes have uncorrelated errors. The A = 
case is of special interest because it reveals which (i.e., 
which frequency/radius pairs) are so constrained by the 
data that even without regularization they may be reliably 
determined. In this limit, ^pj- = '^ij^ y so the uncertainty in 

li is given by Sif = W~^^ where W^j = Y^k TkiTkj / {(Tiaj). 

can be computed by singular value decomposition; 
in practice, we replace the singular values with a small 
number. We show an example of this procedure in Figure 
9. 

The "formal accuracy" of our inversion is illustrated by a 
plot of Ui = II /5Ii (the superscript r stands for recovered 
intensity) as a function of frequency and radius, computed 
from equation (Figure 9). The results for both A = 
and the maximum A consistent with the data are shown. 
In the case A = 0, we use the original h instead of the 
recovered values to compute Ui. The formal accuracy de- 
pends on the true surface brightness: for a given radius, 
Ui tends to peak where the flux is largest. Also, Ui dimin- 
ishes at large radii, since those radii aren't monitored for 
long enough to truly determine li. 

In real solutions A 7^ 0, and the smoothing operator 
correlates the intensities at neighboring radii sharing the 
same frequency. The uncertainties in this case are most 



easily estimated by a Monte Carlo procedure in which the 
lightcurves are perturbed by random realizations of noise 
in the data. The distribution of li after re-solving each of 
these realizations gives the random error in li. In evalu- 
ating these estimates, it is important to understand that 
points weakly constrained by the data have little sensi- 
tivity to measurement errors because they are primarily 
determined by the smoothness constraint. As a result, 
their random errors are artificially small. 

To check for systematic errors in the intensities, we will 
compute the difference between the original and recovered 
surface brightness, A/^ = 1° — II, where superscript o 
stands for "original." The systematic error is, of course, far 
more strongly model-dependent than the random error. It 
depends on the real intensity distribution, the character 
of the data (particularly the sampling), and the inversion 
scheme. T hese considerations will be discussed at greater 
length in § |3l[ 

We estimate the uncertainties in the model parameters 
C, two ways: through the same Monte Carlo procedure as 
for the li, and also through mapping out the x^ found by 
direct solution of the original data (at fixed A) for different 
choices of C,- Examples will be shown in §^ 



3.3 



3. SIMULATIONS 

3.1. Range of Parameters Examined 

To determine how the inversion works in practice, we 
performed simulated inversions, varying the parameters 
describing the underlying model, the parameters describ- 
ing the data set, and the parameters specifying the de- 
tails of the solution technique. By varying the model pa- 
rameters, we learn about whether the method is sensitive 
to the intrinsic nature of the quasar, or the microlensing 
event; by varying the observational parameters, we deter- 
mine what the requirements will be for successful experi- 
ments; by varying the solution parameters, we learn how 
to tune the solution technique for optimum results. 

3.1.1. Model parameters C, 

In all of our simulations we assume that the intrinsic 
radiated intensity in the fluid frame is a black body at 
the local effective temperature, isotropic in the outer half- 
space. Detailed non-LTE spectra computed for our fiducial 
parameters (see the next several paragraphs) are consis- 
tent with both the observed optical/ultraviolet spectrum 
and the microlensing size constraint for the Einstein Cross, 
given the freedom to choose an extinction correction and 
macrolens magnification (Hubeny & Agol, in preparation) . 
However, we examined a number of possibilities for the 
other parameters defining the intrinsic character of the 
quasar and the microlensing events. 

Although we do not know the mass of the black hole, 
we may set reasonable bounds on what it could be. If the 
intrinsic bolometric luminosity is 3 x 10''^ erg/s (Ranch 
& Blandford 1991), the quasar would be at its Eddington 
hmit if Mbh ~ 2 x IQ^Mq. On the other hand, the size 
of the optical emitting region is limited to no more than 
~ 2 X 10^^ cm at ~ 10^^ Hz (quasar rest frame) from mi- 
crolensing (Wambsganss et al. 1990). If this equals the ra- 
dius of maximum emission for an accretion disk (~ lOr^), 
then Mbh ^ 10'^ Mq. On this basis we suppose that the 
true mass is between 2 x lO^M© and lO^M©. Our choice 
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for the fiducial model will be 10^ Mq. 

The units we used in the simulations are Vg for length 
and At, the sampling rate, for time. The units of Vc are 
then rg/At: 

where Vc is the caustic velocity in km/s with distance mea- 
sured at the quasar plane, while time is measured at the 
observer. The caustic velocity, Vc, is quite imcertain, but 
is likely to be in the range 3000 — 5000 km/s (Wyithe et al. 
1999). For At = 3 days, this corresponds to 0.5 < Uc < 4, 
the range we span in our simulations. In the fiducial 
model, we choose Vc = 1- 

We try values of (l)c between and 27r, with n/2 for our 
fiducial model. 

We choose to = for our fiducial model, and also try 
^0 = 5 to see whether this technique works when the cen- 
tral time of the monitoring does not coincide with passage 
of the caustic line across the center of the black hole. 

There are no observational estimates of the inclination 
angle; however, unification arguments for radio-loud AGN 
suggest that quasars are less face-on than blazars, but 
closer to face-on than radio galaxies, so we choose a fidu- 
cial inclination of 30°. We also look at cases with ^ = 
and 6» = 60°. 

Because accretion can spin up black holes, and be- 
cause Kerr holes permit more efficient accretion than 
Schwarzschild black holes, we choose — 0.998 for our 
fiducial model, but also study one example with a = 0. 

The last parameter is the accretion rate. With the fidu- 
cial choices for the other parameters, the observed spec- 
trum is best reproduced with m w 1 = M/{lMQ/yr). 

3.1.2. Observational parameters 

We vary the number of observations, time sampling in- 
terval, SNR, and number of observed wave bands, as well 
as the model parameters. It is especially important to 
determine how the quality of the result depends on the 
number of observations because these observations must 
be targets of opportunity, and thus will impact other ob- 
servations at a given telescope. We explore what happens 
for experiments with between 5 and 41 observations (in 
all cases, we assume uniform spacing). The ratio between 
At and the duration of the microlensing event is implicitly 
given by Vc- 

We try two choices for the SNR (as measured outside the 
microlensing event): 100 for each image in the best case, 
and 50 in the worst case (these were chosen based on cur- 
rent ground-based errors, Rachel Webster, priv. comm.). 

Since disks are broad-band emitters, a broad range of 
observing frequencies is necessary. Observations in the 
four wave bands V, B, R, and I (or equivalent) should be 
routine; observations in U, J, H, K, or in the UV with HST 
will be much more difficult to obtain, but will yield much 
more information. To see just how important the addi- 
tional bands are, we try using just ground-based data in 4 
or 8 bands, or 8 ground -I- 3 HST bands in the best case. 
The short-wavelength bands are especially important for 
hot disks, since the deepest part of the potential well is 
seen at the shortest wavelengths. 



3.1.3. Solution parameters 

Several considerations determine the number of fre- 
quency and radius points at which we may solve for the 
surface brightness. The number of frequency points is not 
simply equal to the number of colors at which the quasar is 
monitored because of the extensive Doppler shifting. We 
found that in practice the best solution grids in both fre- 
quency and radius space were logarithmic. We solve for 
the intensity at frequencies equally spaced logarithmically 
between 3 x 10^^ Hz and 10^^ Hz, and radii equally spaced 
logarithmically from rg to rout, with rout = 500rg. In the 
initial testing of the inversion using the fiducial model, we 
found the smallest number of radii and emitted frequen- 
cies for which we could obtain = Np for some A was 15 
radii and 10 frequencies, which we subsequently used for 
all the simulations. 

3.2. Best-case simulation 

For our "best-case" simulation (designated Al in Ta- 
ble 1), we fixed ^ at the fiducial parameter choices. The 
observational parameters were: 41 observations, SNR = 
100, and 11 spectral bands. 
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time (At) 

Fig. 5. Solid points are lightcurves with noise added. 
The top two curves have been shifted upwards by the 
amount indicated for clarity. The solid lines are the 
lightcurves from the reconstructed disk profile. 

Figure 5 shows the lightcurves for this example, with the 
observed bands de-redshifted for z — 1.695. Note that the 
higher frequencies, which come from nearer to the black 
hole, are magnified more strongly than the lower frequen- 
cies. Figure 6 shows the original and reconstructed disk 
intensity, I{vc,rc), as a function of frequency and radius 
for the best fit parameters for this best case (see Table 
1). The overall shapes arc reproduced quite well. Figure 7 
shows the same results in a different format: we have mul- 
tiplied surface brightness times r^, and plotted the data 
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differently for clarity. Note that the low and high frequen- 
cies and small and large radii are poorly constrained since 
the simulated lightcurve only covers — 20rg < r < 20rg 
and 1600 A < A < 3 ^m. Consequently, at these points 
the regularization tries to make the flux per log radius con- 
stant as a function of radius. Figure 8 shows the ratio of 
the reconstructed to the original one-dimensional profile, 
P^^{x) (computed from I^) for run Al. 



U ,0 
( erg/cm Vs/Hz) 




Fig. 6. Specific intensity I^ire) versus emitted frequency 
and radius. The shaded surface is the original surface 
brightness, while the skeleton plot is the surface bright- 
ness recovered. 




ratio 




Fig. 8. Ratio of recovered to original one dimensional 
disk profile as a function of position and observed fre- 
quency. 

To discuss the reliability of this solution, we begin by 
contrasting the region in the r^-Ve plane where the ran- 
dom error is predicted to be relatively small with the re- 
gion where the actual error is small. As can be seen in 
Figure 9, the region of large Ui for A ^ largely, but 
not entirely, coincides with the region of large /°/|A/i|. 
Moreover, both of these regions follow a track defined, not 
surprisingly, by the requirement that r^I^^ (re) is relatively 
large. Elsewhere in the plane, the contribution to the flux 
is so small that the intensity is virtually unconstrained by 
the data. 




log(ryr ) 



Fig. 9. Plot of /f/Ai (a) and formal accuracy, Ui^ of 
reconstructed disk profile for A — 1.1 (b) and A = (c). 
The dotted line in panel (b) shows where the peak of B^, 
occurs at each radius. 



Fig. 7. Specific intensity I^ive) times versus Ve for 
each emitted frequency. The solid curves are the original 
/i^Tg, while the dashed curves are the minimum and maxi- 
mum values of the recovered I^rl for 20 simulations. Each 
curve is shifted upwards by 2 with respect to the curve be- 
low - the zero point is for the lowest curve. We have not 
plotted negative intensities. 



Where Ui ~ /°/|A/i|, the error is predominantly ran- 
dom error, and 5Ii is a good predictor of its magnitude (in 
fact, in this region A/^ has a Gaussian distribution of the 
correct width). However, there is also a zone on the large 
radius side of the high-intensity track where the systematic 
error is as large or larger than the random error. The na- 
ture of this systematic error is revealed by studying Figure 
7. The smoothing condition tends to raise the intensity in 
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regions where it should be small, and diminish it where it 
is large. Because Ui is rarely large enough to be interest- 
ing where li is small, it is the latter effect that dominates 
in the region of the r-g-z/e plane highlighted in figure 9. 
At least within the context of this model, this systematic 
error is not the result of the specific choice of smooth- 
ing constraint: wc have tried a second-order linearization 
scheme to see if we could get rid of the systematic devia- 
tion; however, we still found that the recovered intensity 
was flatter than the original. If we relax the condition that 
the intensity should be zero at the last radial bin, then the 
intensity approaches a constant for each frequency at large 
radius. 



ratio 




i/(Hz) 



Fig. 10. Ratio of recovered to original one dimensional 
disk profile as a function of position and observed fre- 
quency, using the technique of Grieger et al. (1990). 

We performed a regularized inversion using the Grieger 
et al. (1991) technique for comparison. Figure 10 shows 
their inversion on a data set equivalent to run Al. Since 
the regularization constraint attempts to smooth P^, the 
Doppler peaks are smoothed over, and the noise from the 
lightcurve appears to still be present in the P^. 

After maximizing A consistent with — ^f, we 
then fix A and vary each component of Ci minimizing 
with respect to the other parameters. This proce- 
dure shows how well each parameter can be constrained 
for a given simulation, or whether there are other local 
minima. In Figure 11 we show the Ax^ for each model 
parameter. For this particular model, the physically in- 
teresting parameters, inclination angle (/i) and caustic 
velocity (wc), have well-defined minima. The time of 
origin crossing (io) and the caustic crossing angle (0c) 
are also well-behaved. The black hole spin has a rather 
flat distribution. However, the minimum does lie 

at the correct value. Figure 11 also has a histogram of 
the parameters from each noise realization, showing that 
the minimum of the x^ distribution has few outliers. 
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Fig. 11. Change in x^ versus each parameter for the Al 
case. The vertical lines show the values of the original 
parameters. The histograms show the best-flt parameter 
results of 20 lightcurve realizations, with the right hand 
axis labelling the number in each bin. 

3.3. Varying model parameters 

Table 1 shows the results of varying the model param- 
eters, keeping the observational parameters fixed at the 
"best-case" values. 

Runs Al-AlO each have recovered surface brightness 
which look similar to Al, and the Ui are quite similar. 
Runs Alia and A12 reproduce the intensities well for the 
lowest frequencies and for radii outside rms] however, the 
recovered intensities are non-zero inside r^s, contrary to 
the input model. 

Table 1 shows the average and standard deviation of the 
recovered parameters, C, measured for 20 Monte Carlo re- 
alizations. In all the cases we examined, the distribution 
of recovered parameters is centered near the true model 
parameters, showing that there are no systematic offsets 
introduced by our inversion. This is encouraging since it 
means that this technique has the potential to measure im- 
portant global properties of the accretion disk/black hole 
system. 

Run A2 shows that we can determine the time that the 
caustic crosses the black hole rather accurately. We have 
also tried cases with — ±10, and we find that these are 
also measured quite well, and that the intensity is repro- 
duced as well as in the Al run. 

Runs A3 and A4 show that we can distinguish between 
different caustic velocities, which means that we can con- 
strain the black hole mass in terms of the lens velocity. 
Runs A5 and A6 show that we can measure the disk incli- 
nation angle for a wide range of intrinsic angles. Runs A7, 
A8, A9, and AlO show that we can measure the angle at 
which the caustic crosses the accretion disk rather accu- 
rately. In some cases there is a degeneracy between 0c and 
27r — 4>c when the disk inclination is small, but this should 
not affect the recovered surface brightness since the disk 
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is approximately symmetric in this case. 

Figure 12 shows the topology for each parameter for 
Run AlO, a somewhat special case in which <f>c = tt. The 
parameters Iq and have local minima away from the cor- 
rect minimum; however, these can be ruled out because 
some inferred intensities have large negative excursions in 
the false minimum. 
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Fig. 12. Plot of Ax^ (minimized over all other parame- 
ters) vs. each parameter for run AlO. The solid vertical 
lines show the original parameters. The dotted lines show 
the absolute value of the sum of the inverted intensities 
which are negative (the scale is from to 1). 

Constraining the spin can be difficult, particularly when 
0c — For example, in Run AlO Og is not constrained at 
all (see Figure 12). has as deep a minimum at = as 
it does at the correct value as = 0.998. Only if the coarse 
search in C, space is lucky enough to discover the true mini- 
mum will the Levenberg-Marquardt procedure home in on 
the correct value. It is not clear why some 4>c are more 
favorable for determining Og. 

Whether Os can be constrained at all depends on how 
strongly one believes in the model. If no emission is per- 
mitted inside the marginally stable orbit (Run Allb), Os 
can be constrained because there is a sizable difference be- 
tween the marginally stable orbit around a Schwarzschild 
black hole (6rg) and a maximal Kerr black hole (~ Vg). 
However, if one is unwilling to make this assumption, 
the distinction between the spins largely disappears (Runs 
Alia and A12). The reason for this indistinguishability is 
shown in Figure 13, which shows a contour plot for the 
redshift as a function of position for black holes with spins 
Os — 0.01 and = 0.99, including the regions inside Vms- 
The two plots are almost identical around lOrg, where 
most of the observed radiation comes from in this model. 
Figure 14 shows the x^ (minimized over all other param- 
eters) vs. Os and recovered paramters for 25 simulations 



assuming that emission only occurs outside of Tms- The 
X^ has a clear minimum near the correct spin, and the 
simulations show that the spin can be rather accurately 
recovered. Since Vms increases with decreasing spin, we 
can only hope to obtain a lower limit on Os by assuming 
f > ^ms- Indeed, the x^ vs. Os is flat in the case of zero 
spin (Alia). 




Fig. 13. Redshift as a function of position for 
black holes with spin = 0.01 (solid line) 

and fls — 0.99 (dashed line). 




Fig. 14. Plot of the x^ vs. as (minimized over all other 
parameters) for cases Al and AlO, but with the additional 
assumption that emission only occurs outside Tms- The 
solid dots show x^ for 0c = 7''/2 (case Al) while the open 
circles for AlO. The solid histogram shows the resulting 
Os measured for 25 monte carlo simulations for Al; the 
dotted line for AlO. 

If the disk is much hotter than we have assumed, then 
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the parameters will not be as well constrained as those that 
we have used, as we would then sample only the outer re- 
gions of the disk. To illustrate this point, we have run a 
somewhat unrealistic model, A13, with the fiducial accre- 
tion rate, but a black hole mass of 2 x 10®Mq (near the 
Eddington limit). A standard blackbody accretion disk 
around a black hole with this mass cannot fit the obser- 
vations as its spectrum is too steep and the magnification 
must be much larger than in standard models of the lens 
galaxy. The error on the measured spin is much larger 
than for Run Al (see Table 2). The Ui for Run A13 are 
comparable to those in Run Al, so the intensities are re- 
covered similarly well. 

To see how well we can perform the inversion when the 
assumption of smooth radial variation is incorrect, we mul- 
tiplied the accretion disk intensity by 
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log(r) - log(r„) 



log(r'o„t) - log(r„) 



(13) 



which makes the disk three logarithmically spaced annuli 
(runs A14 and A15). Surprisingly, the recovered model 
parameters, C are accurate. Whether the radial varia- 
tions can be discovered depends on the number of obser- 
vations. In Run A14 (15 observations), the correct overall 
shape is found, but the radial modulation not reproduced; 
in Rim A15 (41 observations), the radial dependence of 
the recovered intensities is more nearly correct. This indi- 
cates that the inversion is only accurate if our smoothness 
model assumption is met on the smallest scale probed by 
the sampling. Figure 15 shows the results of run A14. 




log(ryr|,) 

Fig. 15. Plot of original (solid lines) and recovered 
(dotted lines) for blackbody disk with fluctuations added 
(model A14). Each curve shows a different frequency (from 
3 X 10"'^'* to 10^^ Hz), shifted by 2 units for clear separa- 
tion. The dotted lines show the maximum and minimum 
recovered intensity from 20 Monte Carlo realizations, with 
only positive intensities plotted. 



We have not tried breaking the assumption of azimuthal 
symmetry, as the transfer function is computed assuming 
it. 



3.4. Varying observation parameters 

The question of how many observations are necessary is 
addressed with Runs M2a-d (Table 2). First, we compare 
fewer observations (21) at the same sampling rate (M2a) 
and for the same duration (M2b). In each case, the rms 
scatter of the model parameters is remarkably small com- 
pared to Run Al. Thus, if the underlying model is correct, 
it appears that we can determine the model parameters 
with a high degree of accuracy with relatively few obser- 
vations. 

A smaller number of observations, however, impairs our 
ability to reconstruct the true surface brightness of the ac- 
cretion disk, as the inversion relies more upon the smooth- 
ing constraint than on the actual data. This can be seen 
in Figure 16, which shows the derived formal accuracies 
for various numbers of observations. 
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Fig. 16. This shows a plot of the Ui as a function of 
radius for various numbers of observations at a single fre- 
quency: solid line is Nt = 41 (Al), dotted line is Nt = 21 
(M2b), short-dashed line is iVt = 11 (M2d), and long- 
dashed line is iVt = 5 (M2e). 

The total amount of information about the intensities 
that may be gleaned also depends on the number of ob- 
servations. It is obvious that the number of frequencies at 
which the intensity may be inferred scales in proportion 
to the number of bands whose lightcurves are measured. 
In addition, comparing Al with M2b-d shows that the 
number of radial points with reliable solutions increases 
slowly with increasing number of observation times Nt. 
This is because each observation constrains most strongly 
the minimum radius where the caustic crosses. Since the 
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observations are spaced linearly, while the radii are spaced 
logarithmically, the number of radial points constrained 
by the data is oc log{rmax / rmin) oc log(T/Af) oc log(iVt), 
where T is the total duration and rmax and r„j„ are radii 
corresponding to the radial limits of the region with a re- 

1/2 

liable solution. In addition, Ui oc N^' (see Figure 16), in 
the usual fashion of signal-to-noise ratios. 

In addition to the number of observations, the SNR of 
each observation afi'ects the quality of the inversion. Runs 
M3a-b with SNR=50 demonstrate this dependence. For 
Run M3a, the recovered parameters have similar errors to 
those in Run Al, except for the errors on the spin. For 
Run M3b, the errors on recovered parameters arc roughly 
double those in Run M2d. For fixed model parameters 
and fixed number of observations, on average the errors 
on I are directly proportional to the errors on F; in other 
words, (Ui) oc SNR, where the average is over 20 simu- 
lations. For individual simulations, the Ui depend on A, 
which causes scatter in Ui. 

Next, we look at the dependence of the solution on which 
frequencies are observed. Run M4a shows that the best fit 
model parameters, C are not strongly dependent on the 
infrared frequencies. However, excluding infrared bands 
reduces Ui at the lowest frequencies. Run M4b shows that 
the best fit ^ are strongly dependent on observing the ul- 
traviolet frequencies. This run has the largest errors on 
model parameters and poorest agreement with the true 
parameters. In addition, the Ui at the highest frequencies 
are strongly reduced for this run. Comparing Runs Al, 
M4a,b, we find that the number of frequencies (at the ac- 
cretion disk) constrained by the data oc Ni,, which shows 
why broad frequency coverage is crucial for mapping out 
the disk spectrum. 

Runs Al, A3, and A4 show that both the number of 
radial points with reliable recovered intensities and Ui are 
insensitive to the total duration of the monitoring within 
the range of uncertainties in Vc- Runs A2 shows that the 
time at which we begin the observations does not strongly 
affect the derived parameters or surface intensity, as long 
as we catch the region near the peak. 

4. DISCUSSION 
4.1. Observational requirements 

Experimental design depends on which questions are the 
goal, and to what accuracy one aims to answer them. Con- 
sequently, this discussion, much like our earlier discussion 
of the experimental errors, divides according to whether 
the primary aim is to map the surface brightness of the 
disk, or to infer the model parameters. 

If the goal is to obtain a map of the the disk intensity as a 
function of radius and frequency over as large region in the 
Tg-i^e plane as is possible, it pays most to invest in multiple 
monitoring bands because the number of radius-frequency 
pairs for which it is possible to find a solution is oc , but 
rises only logarithmically with Nt. We expect the intensity 
to vary as a power law in radius, so a better observation 
strategy might be to space observations logarithmically in 
time; however, this will be difficult to achieve in practice, 
which is why we have assumed equal time spacing. The ac- 
curacy of the solution is directly proportional to the SNR 
in the data. To map the regions nearest to the black hole, 
the highest UV frequencies are crucial (although this de- 



pends on the assumption that the spectral peak moves to 
higher frequencies at snialk;r radii): however, lower fre- 
quencies are required if one is interested in obtaining a 
broad-band spectrum. Note also that a higher sampling 
frequency increases the chance that one will obtain moni- 
toring data during the time when the caustic line is near 
the black hole. To obtain high formal accuracy for a broad 
range in radius requires a lightcurve that is finely sampled 
for a long time period. To determine the intensities, the 
model parameters must also be well-defined. When the 
observations are too few or the SNR too small, the uncer- 
tainties in ^ contribute to the uncertainty in Jj. 

If the goal is to simply constrain the model parameters, 
and one believes that the model is correct, then only a 
small number of observations may be required (5 obser- 
vations were sufficient in our simulations to constrain all 
model parameters to better than 20% at la; 11 observa- 
tions for < 10% la accuracy). Since the black hole mass 
and caustic velocity are quite uncertain, it may be neces- 
sary to have more observations to be sure to obtain the 
few near the peak of the caustic crossing that are essential 
for success. The highest frequencies are most important 
for constraining model parameters. 

In our simulations we assumed that only statistical er- 
rors affected the fluxes, and that these errors can be es- 
timated from the observational data. There might also 
be systematic errors in the fluxes due, for example, to 
emission line contributions, inaccurate calibration between 
bands, or cosmic rays. The quality of the result could also 
be affected by these problems. 

4.2. Under what circumstances might the model fail? 

The inversion scheme we propose assumes a model that 
is plausible — a geometrically thin relativistic accretion 
disk in which azimuthal variations are quickly smoothed 
out, microlensed by a caustic system whose basic length 
scale is much larger than the size of the bright region of 
the disk. However, we are by no means guaranteed that 
even this general framework is correct. In this section we 
discuss what would happen if our method is attempted, 
but one of these assumptions is invalid. 

The most basic of our assumptions is that the surface 
brightness varies smoothly as a function of radius. Given 
sufficiently dense sampling with good SNR, even quite 
sharp gradients could be recognized by our procedure. On 
the other hand, if the sampling is inadequate, the existence 
of such features would appear only as a troubling inability 
to find a solution with adequate x^; unless A is taken to 
be very small. 

Another potential source of trouble is departures from 
azimuthal symmetry. In the absence of microlensing, 
"spots" can modulate the lightcurve on the orbital pe- 
riod if they are in the relativistic portion of the disk 
and the inchnation is relatively large (Abramowicz et 
al. 1991). In this case, that would mean periods of ~ 
31Af9(r/10rg)3/2(l + z)/2.7 days (in the observed frame). 
Because the largest observed variations in the Einstein 
Cross are ~ 10% on year-long timescales, this effect cannot 
be too strong in this system. However, there might be a 
range of spot brightness in which they are too weak to show 
up in the ordinary lightcurve, yet strong enough to cause 
some periodic modulation of the lightcurve during a mi- 
crolensing event (Gould & Miralda-Escude 1997). Because 
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the likely duration of a microlensing event (a few weeks) 
is several to ten tirncis the orbital period for the brightest 
part of the disk, it is possible that this effect might be seen 
directly in the lightcurve. If not, they might still make it 
difficult to find a solution with acceptable x^. Particu- 
larly if the break in azimuthal symmetry is approximately 
(X e"" and the disk is nearly face-on, a strong spot might 
be confused with Doppler boosting, leading to a mistaken 
inference for the disk inclination. 

Several physical effects might make disks geometrically 
thick — radiation pressure support if the luminosity ap- 
proaches Eddington (Abramowicz et al. 1988), gas pres- 
sure support if the ions retain most of their heat (Rees et 
al. 1982, Narayan & Yi 1995), or an optically thick out- 
flow. If any of these mechanisms acts, the orbital velocity 
at the photosphere would no longer be that correspond- 
ing to circular free-fall in the equatorial plane of the black 
hole, so that the general relativistic transfer function we 
apply would no longer be valid. We would expect, then, 
difficulties in finding a solution with acceptable . but 
the portion of our solution describing the outer regions of 
the disk should be only weakly affected. 

The surface brightness model we have used for the in- 
version simulations is roughly consistent with both the 
observed spectrum (given the uncertainty in reddenning) 
and the current microlensing size constraint, but it is not 
unique. If the actual spectrum has strong emission in the 
far ultraviolet (which could be true if the reddening is 
greater than our estimate), the observable portion of the 
spectrum will be dominated by emission far from the hori- 
zon, weakening all the relativistic effects. If so, the disk 
inclination and black hole spin will be poorly constrained. 

Uncertain reddening can have other effects, also. Be- 
cause we can expect it to be uniform across the face of 
the disk, it should not affect the inferred radial profile 
of the disk, but it could well introduce additional uncer- 
tainty into our estimate of the intrinsic disk spectrum at 
any given radius. 

Our assumption that the bright part of the disk is small 
relative to the caustic length scale is unlikely to be bro- 
ken, except in the outer regions of the accretion disk. In 
this case, for which the microlensing optical depth is ^ 1 
(Witt & Mao 1994), the caustic scale is essentially the size 
of the Einstein ring due to a single star. Consequently, 
the ratio between the disk size and the caustic scale is 
only ~ 0.01Af9(r/10rg)m~i/2(/i/0.75)i/2^ where m is the 
mean mass of microlensing stars in Solar units, so accord- 
ing to the Grieger et al. (1988) criterion, the caustic as- 
sumption is likely to be valid out to 100 Vg . In our fiducial 
model, for example, the flux at lOOr^ peaks at an observed 
wavelength 1[(1 + z)/2.7]/im, where we have scaled to the 
redshift of the Einstein Cross. 

However, there might be difficulties in practice from a 
related problem: measuring AqF^. Again, if the disk size 
is much smaller than the Einstein radius of a single star, 
then ^0 if' approximately constant during a high amplifi- 
cation event; thus, the same criterion for success applies as 
in the previous paragraph. Of course, observations when 
the quasar is outside the caustic are still required to mea- 
sure Aq. Given the expectation that the quasar is smaller 
at higher frequencies, observing at the highest frequencies 
will provide the best constraint on Aq. To determine how 
these difficulties and those listed in §1.2 will affect the 



inversion, we are currently running simulations of full mi- 
crolenscd lightcurves appropriate for the Einstein Cross, 
to which we will apply our inversion algorithm (Wyithe & 
Agol, in preparation). 

In any of these instances of model inappropriateness, the 
impact on specific inferred parameters depends somewhat 
on details of the inversion procedure. For example, if the 
procedure we have outlined is followed (i.e., minimizing 
by varying C, at fixed A, then raising A until x^ meets 
our definition of acceptability), difficulty in achieving sat- 
isfactory x^ reduces the ultimate A. This means that, in 
effect, more degrees of freedom are "spent" on fitting the 
J(r, v), leaving fewer for defining C,. If the problem is lack 
of smoothness in the radial profile, this transfer of effort is 
reasonable; if the problem is different, however, and if one 
cares about the accuracy of the ^ parameters, one might 
choose to modify the procedure in a way that keeps A fixed 
at a relatively large value. 

One disadvantage of using linear regularization is that 
the intensities are not required to be positive definite, 
though it is impossible to emit a negative number of pho- 
tons. We have tried to incorporate this by trying three 
other methods: maximum entropy; replacing li by log(Ii) 
in B; and the method of projections onto convex sets. Each 
technique finds solutions that are local minima with large 
X^. Another useful technique might be to make the further 
assumption that the spectrum at each radius can be de- 
scribed as a blackbody, and then solve for the temperature 
as a function of radius. These methods are all non-linear, 
and thus intrinsically slow. They therefore impede the ex- 
ploration of parameter space, but might be interesting for 
future work. 

4.3. Multiple microlensing events 

Of the five parameters in ^, three {to, Vc, and (j)c) will 
change from one event to the next, but the other two (Og 
and i) should remain fixed. Since we expect roughly one 
event per year, it should be possible to combine observa- 
tions of several events in order to more tightly constrain 
as and i. 

AA. Connection to X-ray microlensing events 

The optical and ultraviolet continua of quasars are not 
the only portions of the spectrum radiated by the inner 
part of the accretion disk. The X-ray continuum must also 
come from somewhere near that region. It, too, should 
therefore be microlensed in much the same way as the op- 
tical and ultraviolet continuum we have discussed in this 
paper. 

Whether the same technique can be successfully applied 
to the X-ray continuum depends on the same considera- 
tions as discussed in §4.2, but several of them are more 
likely to present problems in the context of X-rays than 
for the optical/ultraviolet continuum. There have been nu- 
merous suggestions, for example, that the X-rays are pro- 
duced in a relatively small number of compact active re- 
gions (e.g., Haardt et al. 1994) that might have substantial 
velocities relative to the disk (Beloborodov 1999). If so, 
the assumptions of azimuthal symmetry, radial smooth- 
ness, and also simple circular orbital motion might all be 
suspect. 

Nonetheless, it would certainly be worthwhile to mon- 
itor the X-ray flux during a microlensing event, in the 
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hope that its emissivity distribution is sufficiently consis- 
tent with our assumptions that it, too, could be mapped. 
Combining this data set with the optical/ultraviolet data 
would also provide an independent constraint on the C pa- 
rameters, which should all be the same for the same event. 

4.5. Summary 

We have demonstrated that monitoring microlensing 
events in the Einstein Cross quasar has great potential 
for both revealing the structure of its continuum emission 
with unprecedented resolution, and potentially constrain- 
ing such basic parameters of the quasar as the spin of its 
black hole (if we assume emission occurs only outside of 
the marginally stable circular orbit), the mass of the black 
hole (modulo the caustic velocity), and the inclination an- 
gle of its disk relative to our line of sight. If this potential is 



realized, and the analytic method we have proposed is im- 
plemented successfully, we may be able to begin answering 
such fundamental questions as: What is the intrinsic local 
spectrum of the disk? How close is it to thermal? And, 
most fundamentally, does the dissipation distribution in 
accretion disks vary with radius in the fashion predicted 
(see equation 1) long-ago? 
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Table 1 

Model parameters for Monte Carlo runs (for = 41, SNR = 100, and N^=11) 





I'O'y'-^'') 


■71 (r / Af) 






ri 

Ug 


\ 

/\ 


Al 





1 


0.86() 


tt/2 


0.998 






-.1±.2 


.!)(i±.()(i 


.84±.()2 


l..")()±.(« 




.01-1.9 


A2 


5 


1 


0.866 


7r/2 


0.998 






-4.8±.2 


.96±.06 


.84±.02 


1.58±.02 


.97±.03 


.002-1.7 


A3 





2 


0.866 


7r/2 


0.998 






-.10±.13 


1.95±.095 


.83±.03 


1.58±.03 


.93±.08 


.01-1.6 


A4 





.5 


0.866 


7r/2 


0.998 






-.18±.2 


.50±.02 


.85±.02 


1.57±.02 


.97±.03 


.001-2.3 


A5 





1 


1.0 


n/a 


0.998 






-.06±.07 


1.02±.04 


.999±.001 


1.53±.08 


.99±.02 


.0001-1.8 


A6 





1 


0.5 


7r/2 


0.998 






-.03±.08 


.95±.05 


.50±.02 


1.57±.01 


.92±.06 


.005-1.8 


A7 





1 


0.866 


7r/4 


0.998 






-.02±.13 


.99±.04 


.84±.01 


.83±.05 


.95±.ll 


.01-1.7 


A8 





1 


0.866 


37r/4 


0.998 






-.03±.l 


1.01±.03 


.86±.01 


2.33±.04 


.99±.02 


.003-1.9 


A9 





1 


0.866 


37r/2 


0.998 






.l±.l 


.98±.06 


.84±.02 


4.72±.02 


.95±.05 


.006-1.9 


AlO 





1 


0.866 


TT 


0.998 






.03±.07 


1.03±.03 


.86±.01 


3.16±.08 


.99±.03 


.005-1.8 


Alia 





1 


0.866 


tt/2 


0. 






-.46±.26 


.95±.07 


.78±.06 


1.57±.02 


.17±.32 


.005-.8 


Allb^ 


-.22±.24 


.87±.07 


.87±.02 


1.57±.02 


.09±.12 


.001-4.4 


A12 





1 


0.866 


7r/2 


0.5 






-.29±.20 


.91±.09 


.83±.07 


1.57±.02 


.55±.25 


.03-.7 



'^We have constrained r > r, 
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Table 2 

Monte Carlo parameters (for to = 0, ^ — .866, (j>c — 7r/2, — .998 model) 



# 


Nt 


SNR 




Vc 


Wc(nieas) 


to{At) 






as 


A 


A13^ 


41 


100 


11 


1 


.93±.07 


-.2±.4 


.85±.04 


1.58±.06 


.89±.19 


.01-2 




15 


100 


11 


2.9 


2.6±.3 


-.04±.06 


.90±.04 


1.57±.03 


.87±.19 


.07-. 9 




41 


100 


11 


1 


.98±.06 


.05±.20 


.88±.05 


1.57±.02 


.79±.20 


.002-. 18 


M2a 


21 


100 


11 


1 


.99±.04 


-.13±.15 


.84±.01 


1.58±.02 


.98±.05 


.5-5.2 


M2b 


21 


100 


11 


2 


1.9±.13 


-.12±.07 


.82±.02 


1.59±.03 


.93±.07 


.2-4.8 


M2c 


15 


100 


11 


2.9 


2.7±.3 


-.09±.l 


.83±.03 


1.59±.03 


.95±.06 


.6-6.6 


M2d 


11 


100 


11 


4 


3.9±.3 


-.03±.07 


.83±.03 


1.56±.03 


.9±.l 


1.7-11 


M2e 


5 


100 


11 


10 


8±2 


0.±.15 


.75±.l 


1.48±.28 


.93±.ll 


9-38 


M3a 


41 


50 


11 


1 


.94±.07 


-.25±.22 


.83±.03 


1.58±.03 


.94±.14 


.006-1.1 


M3b 


11 


50 


11 


4 


3.9±.3 


-.03±.14 


.81±.06 


1.55±.l 


.88±.17 


.8-8.6 


M4a 


15 


100 


8 


2.9 


2.8±.24 


-.04±.08 


.84±.03 


1.57±.03 


.9±.l 


.004-5. 


M4b 


41 


100 


4 


1 


.88±.09 


-.6±.4 


.80±.06 


1.61±.05 


.85±.14 


.25-3.3 



^This run has A/g = 0.2. 

''The intensities have been muhiphed by equation (|l 



